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Abstract. A particular superconducting quantum interference device (SQUID) 
qubit, indicated as double SQUID qubit, can be manipulated by rapidly modifying its 
potential with the application of fast flux pulses. In this system we observe coherent 
oscillations exhibiting non-exponential decay, indicating a non trivial decoherence 
mechanism. Moreover, by tuning the qubit in different conditions (different oscillation 
frequencies) by changing the pulse height, we observe a crossover between two distinct 
decoherence regimes and the existence of an "optimal" point where the qubit is only 
weakly sensitive to intrinsic noise. We find that this behaviour is in agreement 
with a model considering the decoherence caused essentially by low frequency noise 
contributions, and discuss the experimental results and possible issues. 
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1. Introduction 

In the last decade superconducting devices have proved to be promising candidates 
for the implementation of quantum computing [TJ [2]. Single superconducting qubits 
and simple quantum gates have been realized and tested using different schemes and 
solutions with impressive results [3]- [6], providing at the same time a unique framework 
for the study and understanding of intimate aspects of quantum mechanics [7]- [10]. In 
order to further improve the superconducting qubit performances and overcome the 
actual limitations it is fundamental to defeat decoherence due to the variety of solid- 
state noise sources. The first step in this direction is to identify which are the most 
detrimental ones in each specific implementation scheme. This problem has been amply 
investigated in recent years and it is considered as completely understood in the first 
generation of superconding qubits [H] . Those implementations were severely affected by 
low-frequency fluctuations of control variables of different physical origin. Considerable 
improvement has been reached by operating the systems at working points naturally 
protected from low-frequency fluctuations [T2l [13] or by new architectures implementing 
schemes of cavity QED in the solid state, the so called circuit-QED schemes |14j . 
Research along these lines also requires considering innovative materials [T5l EE] • 

In the present article we consider a particular superconducting qubit, the double 
SQUID tunable qubit [I7J HE], where coherent oscillations between flux states are 
obtained by simply applying fast flux pulses (with respect to the typical qubit timescales) 
in the absence of microwaves [19]- [21]. Here we report measurements of coherent 
oscillations at different frequencies obtained by acting on the bias conditions, in 
particular by modifying the pulse height. Interestingly, modifying the oscillation 
frequency also affects its decay time and in particular the shape of the decay envelope. 
With increasing oscillation frequency a crossover between two distinct decoherence 
regimes is observed: an exponential quadratic decay followed by an algebraic behavior. 
In this last regime a working point of minimal decoherence is identified. The existence 
of two decoherence regimes can be explained considering a model system including low- 
frequency fluctuations of the qubit bias and intrinsic parameters. In particular the 
existence of an optimal point of reduced sensitivity to defocusing processes is predicted. 
Our analysis also suggests the possible existence of additional low frequency noise sources 
not included in our model. The paper is organized as follows. In Section 2 the double 
SQUID qubit and its manipulation with flux pulses are described. The experimental 
measurements of the coherent oscillations are reported in Section 3 where we also discuss 
which phase of the manipulation is expected to be more severely influenced by noise 
sources. In the following Section 4 we introduce a theoretical model of decoherence 
processes during the manipulation phase of the experiment and predict the existence 
of optimal points in the present setup. The fit of the experimental data is reported in 
Section 5. We draw our conclusions in the final Section 6. 



Superconducting qubit manipulated by fast pulses 



3 



J l5 Cj h' c 2 



Control ©„ 



Control <!>„ 



<t>. 



(a) 



n 

1 <d' 




Phase 9 



Phase cp 



Phase cp 



Phase cp 



Figure 1. (a) Scheme of the double SQUID, (b-d) Modification of the potential shape 
by changing the bias flux $ c applied to the small loop, (c) Modification of the potential 
symmetry by changing the bias flux <& x applied to the large loop. 



2. The double SQUID qubit manipulated by fast pulses 

The double SQUID qubit consists of a superconducting loop of inductance L interrupted 
by a dc-SQUID, which is a second superconducting loop of smaller inductance I 
interrupted by two Josephson junctions with critical currents and capacitances 
(figure [I] a) . In our case it is L = 85 pH, I = 7 pH, and the junctions are almost identical, 
,h « J 2 = J = 8/iA and C x ~ C 2 = C/2 = 0.4 pF, with {J x - J 2 )/(Ji + J 2 ) ~ 2%. 
The two loops are biased by two distinct magnetic fluxes, indicated as $ x for the large 
loop and <3> c for the small one. When I <C L the system is approximately described by 
a Hamiltonian with a single degree of freedom, the total magnetic flux $ threading the 
large loop or, equivalently, the phase difference across the dc-SQUID (p = (where 
$b = $o/(2vr) = h/(2e) is the reduced flux quantum). The Hamiltonian is the sum of 
a kinetic contribution T = —E c /2d w (E c = (2e) 2 /C is the charging energy, and d w 
stays for the second derivative with respect to the phase), and of the potential 



U{ip) = E L 



(3((p c ) cos ip 



(1) 



where ip x = $ x /®b and ip c = $ c /$b are the reduced bias fluxes, (3(<p c ) = 0o cos(</9 c /2) 
is the tuning parameter controlled by ip c , (3 = 2I L/$ b = Ej/E L and E L = &1/L, 
Ej = 2/o$6 are respectively the inductive and Josephson energy. In our case it is 
E c /h = 0.22 GHz, E L /h = 1920 GHz and Ej/h = 7940 GHz leading to f3 = 4.1. 
Depending on the value of (3(ip c ) , the potential U(ip) has a single well or a double 
well shape, with tp c controlling the barrier height in the double well (figure H]b) and the 
concavity in the single well case (figure[T]c,d), and ip x controlling the potential symmetry 
(figure fib) [221 123] . The described potential presents a periodic behaviour in $ x and $ c 
[19]. In our case the working region is such that $ x ~ and ft (ip c ) is negative. 

In the quasi symmetric condition (\tp x \ 1), fundamental in our scheme, we 
distinguish three important regimes: 

• "W" potential (fig. lb): For /3(tp c ) < —1 the potential U(tp) has double well 
shape (labeled "W" for graphical similarity with this shape), with two minima at 
Vmin § iven b y the solutions of the implicit equation tp^ in + 0(ip c ) sin(^^ in ) = <p x 
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and distant Atp = tp^m — V'mm ■ For fx = (perfect symmetry) the system is 
degenerate and energy levels are arranged in doublets. In practical cases, even a 
small asymmetry \(p x \ <C 1 removes the degeneracy still maintaining the doublet 
structure, with the first two energy levels very close to each other but far away from 
the upper levels (for example, if /3(<p c ) = —1.1 it is sufficient to have |$ z | > 10 _9 $o 
in order to remove the degeneracy). In this case the first two energy eigenstates 
(|0) and |1) ) are flux states in the left and right wells (|L) and \R)). The energy 
gap he = El Atp ip x between them is almost constant in a large range of values of 
(3(cp c ) because of the weak dependence of Atp on (3(<p c ). Approximate analytical 
expressions for the important quantities related to this case have been reported in 



• "V" potential (fig. Id): For 0(ip c ) > —1, U(tp) displays a single well ("V" is used 
for graphical similarity with this shape), and the system can be approximated by 
a harmonic oscillator with level spacing given by 



where Q = 1/v LC. 

• "U" potential (fig.lc): For (3((p c ) ^ -1 we have a rapid transition between the 
double well (W) and the single well (V) case, passing through a quartic potential 
(again "U" is used for graphical similarity with this quartic potential shape). This 
case is particularly interesting for its strong anharmonicity [24] . 

The various shapes of the potential U(tp) and the possibility to easily pass from one 
regime to the other are the fundamental features allowing to use this system as a qubit. 
For example, the qubit "rest state" can be realized when the potential takes the quite 
symmetric double well shape (W). In this regime the magnetic flux states \L) and \R) 
are used as computational states. Qubit manipulations can be performed via flux pulses 
modifying the potential shape, without employing microwave pulses. For example, the 
initial flux state of the qubit can be prepared by strongly unbalancing the potential 
in order to obtain just a single minimum (left or right, according to the unbalancing 
direction), then waiting the time necessary for the complete relaxation in this well, and 
finally returning in the double well condition. A coherent rotation between flux states 
(corresponding to a A rotation in the Bloch sphere) can be realized by applying a fast 
flux pulse on the control tp c in order to change the potential from the double well "W" 
to the single well "V" case ("W-U-V" transition), then waiting in this condition ("V") 
for a time t, and finally returning back ("V-U-W" transition) to the initial double well 
case ("W") where the measurement can be performed. Details of the procedure have 
been reported in refs. [HI 120] . 

Here we summarize the main steps of the manipulation. During the rapid "W-U-V" 
passage, in particular near the "U" Landau-Zener transition occurs equally 

populating the first two energy eigenstates in the "V" harmonic potential, with an 



ref.[23 
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Figure 2. The frequency i7/27r given by cq. © for (p x = as a function of 
(p c = 27r$ c /$o (blue curve). Dashed orizontal lines mark the range of oscillation 
frequencies observed (about 10-20 GHz), corresponding to the range of used top values 
of the applied pulse (defined by the vertical dashed lines). In the inset it is sketched the 
flux pulse used for the qubit manipulation, changing the potential from the two-well 
"W" case to the single- well "V" case (in red). 



initial phase depending on the initial qubit state. Only the first two levels will be 
populated despite of the fact that the final potential "V" is harmonic. In fact, all the 
transitions occur in the strongly anharmonic case "U", where transitions to unwanted 
upper levels can be avoided by accurately choosing the transition speed. During the 
time t spent in the "V" potential, a further phase difference A9 = Qt between the two 
states is acquired. During a final "V-U-W" passage a second Landau-Zener transition, 
specular to the first one, will map the total phase into the probability amplitude for 
the two flux states. For example, if we start the manipulation from an initial left state, 
we will obtain at the end the state = cos(f2£/2)|L) + sm(Qt/2)\R). The rotation 
frequency Q is related to the concavity of the potential in the single well case, and it 
is modulated by the top value of the ip c pulse according to eq. (j2J). The readout of the 
final qubit flux state is performed by a flux discriminator, typically by observing the 
switching of an unshunted dc-SQUID inductively coupled to the large loop of the qubit 



3. Experimental results 

In order to realize experimentally the described operation we need a preparation - 
manipulation - readout sequence. First of all the flux ip x is used to unbalance the 
potential during the preparation phase, and then it is kept fixed very close to zero 
during the remaining time. Then it is applied a flux pulse on tp c presenting a base value 
(pf ase , which maintains the qubit potential in the "W" condition, and a pulse top y^ op , 
which moves the qubit to the "V" condition only for the duration t of the pulse. At 
last the final state is read out by using the dc SQUID discriminator. The preparation - 
manipulation - readout sequence lasts 100/xs in our measurement, and it is repeated N 
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Figure 3. Some experimental oscillations observed for different pulse height. The 
measured frequency is indicated in the top right part of each plot. 

times (with N = 100 -j- 10000, according to the required accuracy) in order to estimate 
the final probability after the manipulation. The sequence is repeated for different pulse 
durations t in order to reconstruct the oscillation of the final probability. Finally, we 
repeat the measurement for different top flux pulse values tp^ op . In order to satisfy the 
condition (3((p^ op ) > —1, the flux yJ op is varied in the interval [0.98 it, 1.12 n]. For the 
sake of clarity we report in figure [2] the dependence of Q on (p^ op . The measurements 
are performed in a dilution refrigerator at 30 mK, on a device based on standard Hypres 
Nb/AlOx/Nb trilayer technology, closed in a system protected by a series of copper, steel 
and mu-metal shields. The dc bias lines are filtered by different LCL low pass filters 
and by thermocoax [26] stages, and attenuators are placed on the fast control lines at 
different temperature stages (see refs. [HI [20] for further details on the manipulation 
and on the setup). 

If the system is prepared in the \L) state, the probability of measuring the qubit in 
the same state after the coherent evolution during a time t in the "V" phase, in the 
ideal case of no noise sources, reads -Pll(£) — [1 + cos(fi£)]/2. In figure [3] we report 
some experimental values of Pll(^) obtained for different pulse heights <£>f op , implying 
different oscillation frequencies (indicated in each panel). Lj] 

Remarkably, changing the pulse height does not only modify the oscillation 
frequency but also influences qualitatively and quantitatively its decay law. This is 
more clearly pointed out considering the envelope of each probability, shown in figure HI 
In contrast with the exponential decay typically originated from quantum noise, the 

\ Similar results have been obtained also by another group in a quite different system [27] . 
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Figure 4. Decay envelopes of the observed oscillations (red points), and relative 
fitting curves (blue continuous lines). 



decay of the probabilities are fitted by the decay law (blue lines in figure H]) | 

m = [i + exp 



(3) 



with independent fitting parameters, 7/ and 777. 

The combination of exponential quadratic and algebraic decay is characteristic of 
defocusing processes due to fluctuactions with 1/ f spectrum of the parameters entering 
the splitting Q [28| 129] . Similar behaviors have in fact been observed in different 
architectures and in the presence of microwave manipulation pulses [29]. The observed 
decay laws therefore suggest that the present experiment is mainly affected by 1/f 
noise in the control fluxes and in the junctions critical currents. Relaxation processes, 
typically leading to exponential decay, seem to be weakly effective. 

The experimentally estimated values of 7/ as a function of the frequency / = 
foy/l + P(ip c ) (/o = f2 /(27r)) are reported in figure (left panel - red dots). We 
observe a regular behavior characterized by a minimum at fo = Qo/(2tt) = 19.8GHz 
[}p T c ° v ~ 7r). This behavior suggests the existence of an optimal point in the present 
setup. The values of 7//, shown in figure [5] (right panel - red dots) are instead quite 
scattered and we observe 777 > 7/, except for frequencies around 10GHz where 777 ~ 7/. 
As a consequence, by changing the operating point, <pJ op , the probability turns from a 
decay approximately exponential to an algebraic behavior for frequencies / close to /o, 
where 7/ is minimum. 

In order to further understand whether this is the appropriate picture, in the 
following Section we theoretically analyze the effect of low and high frequency flux and 
critical current noise in the present setup. Before proceeding with the analysis, we need 
to identify the phases of the manipulation which are presumably more severely affected 
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Figure 5. Red points are the decay rate 7/ (left panel) and 777 (right panel) obtained 
by fitting of the experimental decay curves in figure [3] with eq. ([3]). The blue line in 
the left panel is the fit of the data points with eq. (fT0|) (Section 3). In the right panel 
the blue line is the average value of the scattered values of 7// 



by noise sources. The various parts of the experiment can be addressed separately as 
follows: 

1. Initial phase "W": During the initial "W" phase, when the barrier is very high 
and the system is rigorously prepared in one of the two states, any possible transition 
between wells is blocked by the high barrier. The system remains in an eigenstate to 
high accuracy. 

2. "W-U-V" transition: The effect of noise during this rapid transition is a very 
difficult problem, involving an out-of-equilibrium dynamics. However, for the present 
experiment, noise during this phase is presumably irrelevant. In fact, any noise source 
acting during this transition should affect the oscillations amplitude independently on 
the pulse duration. Therefore we expect that it merely produces a net reduction of the 
visibility of the oscillation at the start time. Since we are able to observe oscillations 
with high visibilities, up to 95% (fig. 2), we conclude that the effect of noise during this 
phase is small. This fact can be tentatively explained by considering that the critical 
region where the change of the potential shape occurs (" U" region) is crossed over in a 
very short time (of the order of picoseconds), shorter than the typical time of coupling 
with the environment (few nanoseconds). 

3. Phase "V": On the contrary, relaxation and decoherence will act heavily during 
the single well condition "V", with an effect depending on the time t of permanence in 
this phase. 

4. "V-U-W" transition: For the return transition from the single to the double well 
potential ("V-U-W") the same considerations of the "W-U-V" transition hold (based 
on the initial visibility of the oscillation), so that one can conclude that also in this case 
the effect of noise is rather small. 

5. Final phase "W" : Finally, when the barrier is raised again in the final condition 
before the readout ("W"), transitions between the two wells are again blocked by the 
barrier. In this case the final state will be a superposition of the left/right states and 
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the effect of dephasing can be quite strong. However, since we are just interested on 
the final population of one of the localized states, the information on the relative phase 
between the two states can be disregarded. 

In conclusion, we expect that the observed decay of the probability can be attributed 
mainly to noise sources acting during the phase "V" , when the system evolves during the 
time t and several repetitions of the protocol are performed. In the following Section we 
will investigate the effect of low-frequency noise during repeated measurement protocols 
in the "V" phase of the present experiment. 

4. Defocusing during repeated measurements 

Superconducting qubits in the various implementations are usually affected by 
broadband and non-monotonic noise |12j |29l [30] . The various noise sources responsible 
for this phenomenology have a qualitative different influence on the system evolution. 
To solve the qubit dissipative dynamics we apply the multi-stage elimination approach 
introduced in Ref. [28]. In simplest cases the effect of noise with large spectral 
components at low frequencies (adiabatic noise), like 1/f noise, and the effect of noise 
acting at frequencies of the order of the qubit splitting (quantum noise), can be treated 
independently. The leading effect of adiabatic noise during repeated measurement 
protocols is defocusing, analogous to inhomogeneous broadening in Nuclear Magnetic 
Resonance (NMR) [31]. It is expressed by the "static path" or static noise approximation 
(SPA) describing the average of signals oscillating at randomly distributed effective 
frequencies [28j EH]- Defocusing can be sensitively reduced by operating at "optimal 
points" where the variance of the stochastic frequency is minimal [32] . Quantum noise is 
instead responsible for relaxation processes. It can be treated by solving a Bloch-Redfield 
master equation, which leads to exponential decay with decoherence time denoted T2 in 
the NMR notation [34] . 

The measurements reported in the previous Section indicate a dominant adiabatic 
noise. Here we concentrate our attention on the effect of noise during the single 
well "V" phase, when the harmonic approximation holds with level spacing given in 
eq.([2]). Defocusing originates from variations of the effective frequency Q during the 
measurement runs which we attribute to fluctuations of the control bias fluxes Stp x , b~ip c 
and of the critical current 5Io, which correspond to fluctuations of (3o, 5(3 . Fluctuations 
of the magnetic fluxes are caused by extrinsic and intrinsic sources. The electromagnetic 
circuit originates quantum noise. In addition, electric lines used to bias the qubit, in 
particular those coming from the room temperature electronics, are responsible for noise 
components mainly at low frequencies, where the filtering is weakly effective (below tens 
of kHz). The intrinsic flux noise acting on the superconducting loops is typically 1/f 
(adiabatic). Different models of flux noise microscopic sources have been proposed, like 
electron hopping between traps either with fixed randomly-oriented spins [35J or with 
spin flips [36]. Spin diffusion along the superconducting surface has also been proposed 
[37] . Fluctuations of the critical current are instead due only to intrinsic causes, like the 
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presence of two state fluctuators in the junction barrier [3HJII5]- In the present article, we 
do not address the important issue of the microscopic source of the fluctuations. Rather, 
we note the existence of low-frequency variations of the effective oscillation frequency 
and describe their effects in the adiabatic approximation. To this end we assume small 
fluctuations of Q given in eq.(|2]) around the nominal values of the control parameters, 
tpT° p and (p x = 0. The expansion of fl to the second order in the fluctuations leads to 
Q + 5Q where 

~ h S(3 + ci 5ip c - \- [ a 2 Sipl + b 2 5/3% + c 2 5ip 2 c + m 2 6ip c 6f3 } . (4) 
i Zq z 

The coefficients of the first and second order terms are (for semplicity here we use (f c 
for <^>) 

1 cos(</? c /2) 1 /3 sin(v9 c /2) 
b\ = = , Ci = = (5) 

+ 4 y/i + PM 

and 

_ PM h cos 2 (y c /2) 

2(1 + /% C ))^ ' ° 2 4(l + /3(^ c ))3/2 



Vl + /3(^) V^ c ' 2(1 + /% c )); ' 2(1 + /%,))^ 

In the SPA the off-diagonal elements of the qubit reduced density matrix (RDM), p(t), 
in the basis of the lowest eigenstates in the "V" potential {|0), |1)} is obtained by 
evaluating the average over the fluctuations 5ip x , 5ip c and 6p 

W(*)> = p 01 (0)e- iQt {e- iSQt ), (8) 

here we assume that fluxes and critical current fluctuations are are uncorrelated random 
variables with Gaussian distribution, zero mean and standard deviations a x , a c and ap 
respectively. We remark that for the fluctuations of the critical current this assumption 
should be checked considering the microscopic source of Iq fluctuations. The power 
spectra of the control flux biases and of the critical current measured in similar flux 
or phase qubits reveal a 1/f behavior at low-frequencies [35j ESj- The variance of the 
corresponding variables, a 2 a [a = x,c,/3q), is proportional to the amplitude of the 1/f 
spectrum, Sa\u) = ira 2 [ln(7 Ma ,/7 ma ) a;] -1 ("y ma and 7m« are the low and the high 
frequency cut-offs of the 1/f region). 

In order to compare with the experiments, we evaluate the average probability (Pll{^))- 
For the chosen initial condition, \ip(t = 0)} = \L) = (|0) + \ l))/\/2 the probability reads 
(-Pll(^)) = \ + (R- e [Poi (£)])• Where we assumed that the sum of the populations of 
the lowest eigenstates of the harmonic potential is constant, poo(^) + pn{t) = 1- This 
assumption is justified in the adiabatic approximation, where populations do not evolve. 
It also holds in the presence of quantum noise, provided leakage from the bi-dimensional 
Hilbert space can be neglected. The averaged probability thus reads 

(P LL (t)} = ^[1 + (cos((JI + 5Q)t))] = ^[1 + d(t) cos(fit + 6)] (9) 
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Figure 6. Variance of Q given by eq. ([T3|) as a function of ip c in the allowed interval 
of values where f3(tp c ) > —1 for a c ~ a x « 3 x 10~ 4 , and different <r^ values 2 x 10 -5 , 
2 x 1CT 4 , 2 x 1CT 3 , 2 x 1CT 2 from bottom to top. 



where d{t) = \(expi5flt)\ and 9 = arg[(exp iSQt)]. The coherence {poi{t)) in the SPA 
is obtained evaluating the average eq.® with 5Q given by eq. fl4]). For short times 
(neglecting contributions oc t A with respect to terms oc t 2 ) d(t) is indeed given by eq.(j3]) 
with 7/ and 777 given by 



li 



Q Jb 2 a 2 



In = ^0 



' /3 + c^c + + rniafa 

Note that the first order terms of the expansion enter 7/, the second order ones 
enter 7//. Finally we include the possibility of incoherent energy exchanges with the 
environment due to quantum noise. In the simplest version of the multi-stage elimination 
approach the overall decay factor of the coherence is 

d{t) = \{expi5nt)\e- t/T2 , (12) 

where the decoherence time due to quantum noise, T 2 , depends on the power spectrum 
of the flux and of the critical current fluctuations at frequency O. 

The existence of an optimal point in the present setup depends on the possibility 
to minimize the standard deviation of the effective splitting Q, E = a/ (5Q 2 ) 
[33] . In fact, from the short-times expansion of eq, 



.0:. 



(10) 
(11) 



(5tt) 2 



([8]) we find \(e- i5nt }\ w 
(St) 2 , thus defocusing is reduced when S is minimal. In the present case, from 



(HI), we find 



2 1 2 
7/ + if 11 



(13) 



Because of the periodicity of the first order terms in (j4j), included in the coefficients 
bi and C\ defined in eq. JS),E is mimimal either at tp c « 0, when a c > crg , or at tp c m tt 
in the opposite case. This is illustrated in figure |6] where we show S for typical values 
of the variances a c and a x j35l [38] and increasing values of ap . Depending on which 
condition applies to a specific setup, the device can be operated at the proper optimal 




Figure 7. Left panel: 7/ (blue) and 7// (red) as a function of ip c for a c ^ <r x ~ 
3 x 10~ 4 , and different ap a values (solid lines for 2 x 10 -5 , dashed lines for 2 x 10 -3 , 
dotted lines for 2 x 10 -2 ). Right panel: 7/ (blue) and 7// (red) as a function of 
/ = CI/2tt in the experimental regime ip c € [0.98 w, 1.2 7r] for the previous values of tr c , 
a x and cr^ = 2 x 10 -5 (solid line), ap = 2 x 10~ 2 (dotted line). 



point 

The dependence of 7/ and 777 on the operating point is reported in figure [7J for the 
same parameters of figure O In the considered regime, 7/ is minimal at the optimal 
point, wherever it is located (left panel). This is a consequence of the fact that under 
these conditions the variance of Q is dominated by the first order terms and £ « 7/. 
In the right panel of figure [7J we show 7/ and 7// as a function of the frequency / when 
ip c varies in the experimental range around <p c ~ 7r (as indicated in figure [2]). In this 
frequency range 7/ is minimal only when ap > a c (corresponding to the dot-dashed and 
dotted curves in figure EJ). We may thus expect that a similar situation occurs in the 
present experiment, with an optimal point at corresponding to a minimum of 

7/. However we note that for the parameters in figure [7J since 7/ > 7//, the coherence 
decay d(t) is dominated by the exponential factor exp(7ft 2 /2), independently of the 
optimal point location. Since in the present experiment we observe a crossover from an 
exponential to an algebraic decay by changing the working point, for some operating 
conditions around tp T c ° v ~ n it should turn out that 7/ < 7//. Since at the considered 
working point (<p x = 0) fluctuations of the flux (p x contribute to the second order in (TJJ, 
the above relation between 77 and 7// may in principle occur when the flux variances 
differ considerably. Such a situation is hovewer unlikely, a dependence of 1/ f flux noise 
on the SQUID size as due to a "global magnetic field noise" has in fact been ruled 
out [35]. Such a situation is also excluded in the present experiment. An illustrative 
case with a x ^> a c is shown in figure Due to the increased weight of the second order 
terms in (jlj), 777 has a local mimimum at the optimal point, (pJ op ~ 7r. However for 
frequencies close to the minimum, being 7/7 < 7/, the decay factor would be exponential 
rather than algebraic, in contrast to the observations reported in the previous Section. 



§ Note that at the optimal point the first derivative of the splitting Q, with respect to either ip c or /?o 
is non-vanishing. It is not possible to reach a optimal point of vanishing differential dispersion as for 
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Figure 8. Left panel: variance E 2 as a function of ip c with the minimum at ps 7r, 
continuous line is (|13[) , the dashed line is the first order approximation E w 7/. Right 
panel: the rates 7/ and 7// as a function of /. The variances in this example have been 
fixed to highlight the regime 77 > 7/7, and are a c w 10 -3 , cj^ ps 10 _1 , crs ps 10~ 4 . 

We conclude from this analysis that the present model, considering adiabatic 
fluctuations of the two fluxes ip c and (p x , and of the critical current, predicts the existence 
of an optimal point of minimal variance of the effective splitting and of minimal 77, 
when first order terms in the expansion of the effective frequency Q dominate. On the 
other side, it is difficult to predict the simultaneous existence of an optimal point at 
^Top n w \y eTe ry x j g minimum and 77 < 777 only close to the minimum. A plausible 
scenario is that the present experiment is sensitive to first order contributions, entering 
77, whereas second order effects entering 777 are masked possibly by additional noise 
sources not included in our analysis. For instance, our model does not include SQUID 
inductance fluctuations with 1/f power spectrum, which are highly correlated to 1/f 
flux noise [39] . 

5. Fit of experimental data: 77, 777 

According to the analysis of the previous Section, whereas for 77 we may expect an 
agreement between eq. ([TO]) and experimental data, we expect that eq.f lTT]) for 777 
has to be supplemented by an additional contribution of different origin. Therefore, in 
order to extract 77, 777 and 1/T 2 from the experimental data we fit the envelopes of the 
experimental curves in figure |3] considering them as independent parameters. 

Indeed, the experimentally estimated 77 for the different oscillation frequencies 
/ can be fitted with eq. OlOp . giving a remarkable agreement with fitting parameters 
o c = 1.4 x 10~ 3 and ap = 0.105, corresponding to er$ c = 223/!$o an d 0% = 1.7 /xA, 
see figure [5] (left panel, blue line). The noise on the flux ip c in addition to the 1/f 
behavior, possibly is also influenced by low frequency noise components due to the 
room temperature instrumentation. For the estimated noise variances, fluctuations 
of the critical current are almost always dominating except for the optimal point at 

charge or flux qubits [IB] . 
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Figure 9. Plot of jj (blue) and jfj = jn + 2.6 ■ 10 9 s 1 (red) as a function of / for 
variances extracted from the fit of 77, <x c « a x « 10~ 3 , crg « 10 _1 . 

/o = fio/(27r) = 19.8GHz where the effect of flux noise emerges. The existence of this 
optimal point is evident from the decays in figure HI where a crossover between two decay 
regimes is observed: a fast decay at lower frequencies is followed by a slower decay at 
higher frequencies. 

On the other side, the value of obtained form the fit is much larger than the 
values reported in other superconducting qubits [531 15H] . An independent check of this 
quantity in the present setup is not possible, since measurements of the low-frequency 
power spectrum of critical current noise are not available at present. One possible 
explanation of the quite high value of a Io is related to the materials used for the qubit 
fabrication. In fact other superconducting qubits based on the same materials and 
fabrication technology and displaying coherence times similar to the ones reported in 
our experiment (few nanoseconds) [21] presented a considerable enhancement of these 
times by improving the used materials, for example by introducing SiN x dielectric films 
instead of standard S1O2 for the crossover wiring [T5J [16]. The full understanding of 
this point requires the repetition of the experiment with the use of different materials 
and technologies. Another possible source of this inconsistency stems from the assumed 
Gaussian distribution for the critical current fluctuations, 5(3 . Evaluating possible 
deviations from the Gaussian approximation requires considering a microscopic model 
of critical current fluctuations. This is another possible extension of the present analysis. 

The values of 777 extracted from the fit of the oscillations in figure [3] are scattered 
around an average value ~ 2.6 • 10 9 s _1 , as shown in figure [5j Equation (llip for 7// 
with o~ x ~ cr c = 1.4 x 10~ 3 and crg = 0.105 predicts 7// <C 7/, in contrast with the 
observations. However, including a constant (frequency independent) noise contribution 
in the second order terms, i.e. defining an effective jfj ~ 7// + 2.6 • 10 9 s _1 , the 
crossover from exponential to algebraic decay where 7/ <C 7// visible in figure [5J can 
be quantitatively reconstructed, as illustrated in figure [9j 

Finally, the reported measurements do not allow for a reliable estimate of the effect 
of quantum noise, included in the exponential decay term with T2. The observations 
are compatible with a decoherence time T 2 with a lower limit of tens of nanoseconds, 
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which is reasonable in this type of qubits [To] . 
6. Conclusions 

In conclusion, we considered the manipulation of a double SQUID qubit by fast flux 
pulses, and observed the decay envelope of the obtained oscillations for different control 
conditions, corresponding to different oscillation frequencies. The shapes of the decay 
envelopes show a peculiar behaviour with a crossover between two distinct regimes. 
These behaviors can be attributed to various sources of adiabatic noise affecting the 
system. The effect of high frequency noise is negligible, and this indicates a correct 
filtering and shielding of the system. 

We observed a crossover between an exponential and an algebraic decay regime, 
with an optimal point where decay is algebraic. We demonstrated that this behavior is 
due to the interplay of first order effects of low-frequency flux and critical current noise. 
In general, intrinsic fluctuation of the critical current dominate except at the optimal 
point where the weaker effect of flux bias fluctuations springs up. The existence of the 
optimal point is an interesting characteristic for possible applications. 

The effect of second order fluctuating terms is still misunderstood and will require 
further investigation. To this end we plan to repeat the experiment with different, 
improved materials. The ensuing noise characterization may possibly provide new 
insights into the low-frequency noise sources in this kind of setup. 
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